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evolution of the shape of superconducting transition R{T) in quasi-one-dimensional 



C I aluminum structures. In nanowires with effective diameter < 15 nm the R{T) de- 

pendences are much wider than predicted by the model of thermally activated phase 
slips. The effect can be explained by quantum fluctuations of the order parameter. 



pH , Negative magnetoresistance is observed in the thinest samples. Experimental results 
O 

O I are in reasonable agreement with existing theoretical models. The effect should have 

Q^' a universal validity indicating a breakdown of zero resistance state in a supercon- 

. ductor below a certain scale. 
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I. INTRODUCTION 



The ability to carry a dissipationless electric current is a fundamental attribute of a super- 
conductor. One might naively expect that this property is preserved at reduced dimensions. 
Unfortunately, inevitable fluctuations result in momentary suppression of superconductiv- 
ity leading to energy dissipation. Being integrated in time, this effect manifests itself as 
finite resistance. Of particular interest are quasi-one-dimensional (quasi-lD) systems, with 
effective diameter a^l"^ smaller than the superconducting coherence length ^, where there 
is only one parallel channel of supercurrent. In such objects, due to fluctuations of the 
order parameter, the effective resistance is never strictly equal to zero. However, a measur- 
able resistance at temperatures well below the critical temperature can be observed only 
in samples with rather small diameters. Early experiments |l| confirmed that in narrow 
superconducting channels the shape of the R{T\ transition can be explained by thermal 
fluctuations Recent experimental [3j - [8j and theoretical [10] - [19[ discoveries claim the 
existence of an additional (non thermal) mechanism of finite resistance in ID systems. The 
issue is of vital importance for the development of superconducting nanoelectronic elements 
designed to carry a dissipationless electric current. 

The finite resistance of a ID channel can be understood that for sufficiently longs system 

of length L ^ ^ there is always the finite probability of a fluctuation to drive instantly a 

fraction of the wire of volume ^cr into a normal state. The energy required for this process 

is the corresponding superconducting condensation energy AF ~ B^^a, where Be is the 

critical magnetic field. If thermal effects solely contribute to the fluctuations, it has been 
ft 

shown [2] that the effective dc voltage V(T, I) is proportional to the probability of these 
events ~ exp(— AF/Zc^T), where ks is the Boltzman constant: 



V{T, I) = n{T, I, L) exp 



AF /2\i/2 J 



sinh ( ^ ) , (1) 



Ic is the temperature-dependent critical current, and Iq = ksT /(po with 0o being the super- 
conducting flux quantum. The current dependence is determined by the sinh{I /2Iq) term 
and the term under the exponent, responsible for the effective reduction of the potential 
barrier AF by the finite current /. An exact form of the pre-factor f2(T, /, L) has been pro- 
posed [2I. Though the function f2(T, /, L) also contains temperature- and current- dependent 
terms, it contributes negligibly compared to the strong exponential dependence. Formally 
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the above process can be described as the thermally activated transition of a superconduct- 
ing system from a local potential minimum to a neighboring one separated by ±27r in phase 
space: so called, thermally activated phase slip (TAPS). One can make a formal analogy 
with a "classical jump" of a particle over the barrier AF stimulated by the thermal energy 
ksT. 

An alternative mechanism has been proposed associated with tunneling through the en- 
ergy barrier AF in phase space jsl- Such tunneling to a state of lower energy, also called 
quantum phase slip (QPS), should provide an additional (to TAPS) channel of energy dis- 
sipation in a current carrying system. For conventional type-I superconductors estimations 
show that measurable deviations from the TAPS mechanism are expected in samples with 
effective diameter 10 nm. There have been a few experiments, claiming observation 

of QPS phenomenon in various superconducting materials: In, Pb and Pb-In tSil, Pb, Sn 

n n n q n 

and Pb-Bi MoGe [5j, Sn and Zn [6[, Al [7j, [8|. Various models consider non-thermal 
mechanisms of the finite resisitivity of a ID superconductor jol - [l9|. In spite of the inten- 
sive research in the field, the matter is far from being settled. From an experimental point 
of view, to some extent, the inhomogeneity of a nanowire might lead to a broad supercon- 
ducting transition R{T) which can be erroneously associated with the QPS mechanism 20|. 
Another source of misinterpretation might come from insufficient filtering from environmen- 
tal RF noise. The resulting parasitic overheating can shift the superconducting transition 
to lower temperatures, distorting the shape of the -R(T) dependence. 

Our motivation was to perform experiments eliminating the uncertainty related to the 
uniqueness of each particular nanostructure. We were able to trace the crossover from the 
TAPS mechanism to QPS in the same sample , successively reducing its diameter between 
measurements, thereby evidencing for a solely size dependent origin of the phenomenon. The 
experiments performed on multiple sets of aluminum nanowires reduced by ion sputtering 
down to sub-10 nm scales showed quantitatively similar behavior. In all structures with 
effective diameter a^/^ < 15 nm the R{T) dependences deviate from the predictions of the 
TAPS model [2|], while the QPS scenario [HI] provides a reasonable fit to the data. 



II. EXPERIMENT 
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In our earlier work 



2l| we demonstrated that low energy Ar~^ ion sputtering progressively 



and non-destructively reduces nanostructure dimensions. The penetration depth of Ar~^ ions 
into an Al matrix at acceleration voltages of ~ 500 eV is about 1.5 nm and is comparable 
to the thickness of naturally formed oxide. The accuracy of the effective diameter cr^/^ 
determination from the normal state resistance, scanning electron microscope (SEM) and 
scanning probe microscope (SPM) measurements is about ± 2 nm. Only those samples which 
showed no obvious geometrical imperfections were used for further experiments. While 
reducing the nanowire diameter, the low energy ion sputtering also provides "polishing" 



removin g su r: 



'ace roughness inevitable after a lift-off process j?! (Fig. 1). In this work the 



method 2l|], 25|] was applied to lithographically fabricated 99.999 % pure Al nanowires with 



initial cross-section of about 100 nm x 100 nm. During measurements the structures were 
immersed into a directly pumped ^He bath with base temperature of about 0.95 K. Contrary 



• 0, B, 



to Refs. m, m, m, we used a four-probe configuration, ehminating necessity to consider 
the contribution of the electrodes. We obtained exact quantitative agreement between data 
measured with dc and ac currents using extensive multistage rf filtering. The results are 
reproducible and do not suffer from hardware-related artifacts. 
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FIG. 1: (a) Atomic force microscope (AFM) image of a typical 5 /xm long aluminum nanowire 
just after lift-off. (b) AFM images showing evolution of the shape of the same nanowire after 

several sessions of ion beam sputtering. Bright color above the horizonatal plane (initial level of 
the substrate) corresponds to metal, dark color below to sputtered Si substrate. Note the reduction 
of the initial surface roughness of the nanowire. 
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III. RESULTS AND DISCUSSION 



After a sequence of sputterings [alternated with R{T) measurements] the wire cross sec- 
tion a is reduced down to < 20 nm, where deviations from TAPS behavior become obvious 
(Fig. 2). Similar results were obtained on several sets of aluminum nanowires with length L 
equal to 1 /im, 5 /im and 10 fim. For larger diameters (e.g., Fig. 3, insert, 55 nm wire) the 
shape of the R(T) dependence can be qualitatively described by the TAPS mechanism. No- 
tice that a size dependent variation of the critical temperature of a nanowire [2^, especially 
pronounced in aluminum results in a broadening of the^T) transition and significantly 
reduces straightforward applicability of the TAPS model 20[. Confidence in assigning our 
results to the manifestation of a new (non-thermal) mechanism comes from the experimen- 



tal observation that in aluminum nanowires (fi 
decrease of the wire diameter (film thickness) 



ms) the critical temperature increases with a 

"1 Q 

22[, [7|. It means that the broadening of the 



R{T) dependencies below the bulk value ~ 1.2 K cannot be explained by a geometrical 
imperfection (e.g., constriction) of an aluminium nanowire. Below a certain limit < 
20 nm fits by the TAPS model fail to provide any reasonable quantitative agreement with 
experiment even assuming the existence of unrealistically narrow constrictions not observed 
by SPM (Fig. 2). 
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FIG. 2: Resistance vs temperature for the thinnest samples obtained by progressive reduction of 
the diameter of the same aluminum nanowire Al-Cul26-3 with length L= 10 /xm. The Langer- 
Ambegaokar-McCumber-Halperin (LAMH) model fitting is shown with dashed lines for 11 and 15 
nm samples with the best fit mean free path i = 3 and 10 nm, correspondingly, Tc = 1.46 K and 
critical magnetic field Bc{0) = 10 mT. Fitting using a simplified short wire model Ref. [11] [Eqs. 2 
and 3] is shown with solid lines. For 11, 12, 13 and 15 nm wires the fitting parameters are: Tc = 1.5 
K; ^ = 0.15, mean free path i = 5.4 nm, 5.8 nm, 7.3 nm, 7.5 nm; and the normal state resitance 
Rn = 7200 kO, 5300 kO, 4200 kO and 2700 kO. 
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Though several theoretical approaches [9j - [19| were proposed to describe the non-TAPS 
mechanism of finite resistance in ID superconducting channels, not all models are suitable 
for direct comparison with experiment, lacking a clear expression for the resistance (voltage) 
as a function of temperature. Where possible, we tried to compare our data with the models 



and found the best agreement with renormalization theory [HI]. The full version of the 
model for systems of an arbitrary length L is rather sophisticated, requiring knowledge of 
parameters not easily deductible from experiment. However, if the wire is short enough that 
only one phase slip event can happen at a time, one can neglect the effects of the interaction 
between the phase slips. In this limit the simplified expression for the rate of QPS activation 
is 

Tqps = ^^^^exp^-Sqps), (2) 

where the action Sqps = A{Rq / ^) / {R^ / L) , with A being a numerical constant, Rq = 
/i/(4e^) = 6.47 kQ and tq ~ /i/A is the characteristic response time of a superconducting 



system which roughly determines the duration of each QPS HI]. The effective (time av- 
eraged) voltage Veff due to the fluctuations can be found using the Josephson relation j^. 
Finally, for the effective resistance of a quasi- ID superconducting nanowire one gets 

Rqps{T) = Kff/I = (3) 

Contrary to TAPS j^, the QPS contribution has a rather weak temperature dependence 
far from the critical temperature and should produce a flnite resistance even at T -^^ for 
sufficiently narrow wires [ll]. One can achieve reasonable agreement between the model 
[Eqs. 2 and 3] and the experiment (Fig. 2). There are four fitting parameters: critical 
temperature Tc, normal state resistance Rn, mean free path i [to re-calculate the dirty limit 
coherence length ^ = 0.85(^o^)^'^^ ]^ and the numerical parameter A of the order of unit 



ll| . The critical temperature and the normal state resistance can be trivially deduced from 
experimental R{T) dependences. Roughly the mean free path i can be estimated from the 
normal state resistivity p, as the product pi = 5 x lO'^^Qm'^ is a well-tabulated value for 
dirty-limit aluminum. In our ultra-narrow nanowires the cross-section is known from SPM 
measurements with ± 2 nm accuracy. Hence, there is some freedom in the selection of the 
mean free path. As a rule of thumb for all our nanowires with effective diameter a^/^ < 20 
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nm the best-fitted mean free patli (at low temperatures) was found to be roughly equal to 
one-half of the diameter (Fig. 2, caption). The observation is quite reasonable, taking into 
consideration that at these scales and temperatures electron scattering is mainly determined 
by the sample's physical boundaries. For all aluminum nanowires the best-fitted value for 
the parameter A was found to be equal to 0.15. As the value of A cannot be calculated within 
the model llj with the required accuracy, we believe that the correspondence between the 



simplified "short wire" model and the experiment can be considered as good. 

No Coulomb blockade has been observed on V{I) characteristics, indicating good ho- 
mogeneity of the wires (absence of tunnel barriers). Below the critical temperature V{I) 
deperidences show ~ sinh{I / Iq) behavior similar to both TAPS [Eq. 1[] and renormaliza- 
tion [ll] models. At sufficiently small currents / ^ 2Jo ~ 20 nA one might expect a linear 
response, while at much higher currents trivial overheating is observed. However, there exist 
an intermediate regime between linear and strongly non-linear regimes: the top part of the 
R{T) transition is not shifted, while the slope of the bottom part decreases with an increase 
of the measuring current (Fig. 3). We associate this peculiarity with the reduction of the 

)otential barrier AF by a finite current / [Eq. 1] , which is qualitatively similar for thermal 

2j and quantum llj mechanisms. 
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FIG. 3: Resistance vs temperature for a a^l'^ = 11 nm it 2 nm sample of aluminum nanowire 
Al-Cull5 with length L= 10 [im measured at various DC and AC currents. Inset: resistance 
vs temperature for three samples obtained by progressive reduction of the diameter of the same 
nanowire. 
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At certain limits the renormalization model [llj predicts the same functional dependence 
of the effective resistance on temperature and current. In the high temperature limit T ^ 
(pol/ks , R ~ T^'*'^^, and in the high current limit / ^ Jq , -R ~ . The dimensionless 
conductance 7 = Rq/Rqp is related to the effective "quasiparticle" resistance i?qp, being 
associated with dissipation provided by the quasiparticle channel. Rqp can be considered 
as a fitting parameter and should be of the order of the normal state resitance Rjq. In the 
metallic phase, when quantum fluctuations dominate the behavior of the system, the residual 
resistance at T — > should correspond to Rqp lli. The results of fitting of our R{T) and 
V{I) dependences to the above relations are presented in Fig. 4. One may obtain satisfactory 
agreement with experimental data, allowing the quasiparticle resistance Rqp to be smaller 
than the normal state resistance Rn of the nanowire. For each sample, the deviation of the 
fitting parameter Rqp from R]s[ is less pronounced for V{I) characteristics. Probably, this 
discrepancy is the result of the experimental difficulty to satisfy the applicability condition 
T ^ (poI/kB, while the temperature should be noticeably below the critical temperature Tc, 
where the whole concept of phase slippage is valid. 
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FIG. 4: (a) Resistance vs. temperature for three samples of the same wire as in Fig. 2 with 
progressively reduced cross-sections. Solid dots are fittings to power dependence R ^ T^i-^, (b) 
V(I) dependences of the same samples taken at close temperatures stabilized with accuracy it 0.1 
mK. Solid lines correspond to proportionality R{T) = V{T)/I ~ In both figures the only 

fitting parameter is the quasiparticle resistance Rqp. 
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Application of a magnetic field perpendicular to the plane of the structures reveals a 
rather unusual effect. At sufficiently small temperatures a pronounced negative magnetore- 
sistence (nMR) is observed (Fig. 5, inset), resulting in a "sharpening" of the R{T) transitions 
in magnetic fields < 25 mT (Fig. 5). nMR is observed only in the thinest samples cr^^^ < 
20 nm and only at temperatures well below the onset of superconductivity. A similar effect 
has been reported for ultrathin ID lead strips {4]. Aluminum is known to be immune to the 
creation of localized magnetic moments provided by the majority of magnetic elements. Our 
samples were e-beam evaporated from 99.999 % pure aluminum target in an UHV chamber 
where magnetic materials have never been processed. Hence, the residual concentration of 
magnetic ions should be negligible. Even if we were to assume their presence, magnetic fields 
where nMR is observed are too small to polarize them. Additionally, we do not see an en- 
hancement of Tc related to the onset of superconductivity, which would be the consequence 
of suppression of the Kondo mechanism by a magnetic field. As the presence of localized 
magnetic moments in our samples is not obvious, a recently developed model of nMR Re 



231] is not applicable. nMR has been predicted in disordered superconducting wires 19 ] 



However, the main result of that model is an exponentially small addition to the TAPS ef- 
fective resistance. This conclusion contradicts our main observation (Fig. 2) and, hence, the 
applicability of the whole model |l^ to our experiments is questionable. For the moment we 
do not have a solid explanation for the nMR effect. The phenomenon might be related to an 
interplay between two field-dependent contributions: AF barrier reduction and suppression 
of the superconducting energy gap A. The first process provides an expected increase of 
the observed resistance, while the second one leads to a reduced quasiparticle resistance Rgp 
due to thermal activation of extra quasiparticles resulting in lower effective resistance R{T). 
There might be a region of magnetic fields where the second mechanism "wins" reducing the 
influence of quantum fluctuations leading to nMR ll(]. Another explanation of the nMR 
might be related to the formation of a charge imbalance region accompanying each phase 
slip event 24] . This non-equilibrium region provides dissipation outside the core of a phase 
slip. The corresponding Ohmic contribution can be effectively suppressed by the magnetic 
field, resulting in nMR. However, so far the validity of the charge imbalance concept was 
only demonstrated at temperatures sufficiently close to and its applicability to QPS at 
temperatures well below the critical one is not obvious. A quantitative comparison with the 
experiments requires further elaboration of the theory. 
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FIG. 5: Slowly recorded (~ 1 h) resistance vs temperature dependence for the 11 nm sample from 
Fig. 2. While sweeping the temperature, a few times the perpendicular magnetic field B = 19.6 
mT was switched on and off. The top branch corresponds to zero field, while the lower one to field 
"on". Inset: resistance vs perpendicular magnetic field for the same sample measured at constant 
temperature and small ac current. 
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IV. CONCLUSIONS 



The homogeneity of the wires is the central point in an interpretation of data related 
to phase slip mechanisms: thermal or quantum. The existence of trivial structural or geo- 
metrical imperfections as constrictions, boundaries and contact regions might broaden the 
i?(T) dependences [20| and can be erroneously interpreted as a "new physic'". The specially 
developed method of ion beam sputtering [2l| , [25|] to a large extent allowed us to study the 
evolution of the size phenomenon, eliminating artifacts related to the uniqueness of samples 
fabricated in independent processing runs. The ion beam treatment polishes the surface 
of the samples, removing the inevitable roughness just after fabrication (Fig. 1). If there 
were no detectable geometrical imperfections in the original (thick) wires, they cannot be 
introduced while reducing the diameter by the low energy ion sputtering. Formally, it can- 
not be excluded that in the original structures there were "hidden" undetectable structural 
defects (e.g., highly resistive grain boundaries) which did not contribute to conductivity, 
being shunted by the bulk. When the diameter of the wire is reduced below a certain scale, 
this type of imperfections might determine the behavior of the quasi-lD system. However, 
the absence of Coulomb blockade and/or gap structure on J — characteristics rules out 
the presence of the highly resistive tunnel barriers. 

The extensive SPM and SEM analysis indicates that after several sessions of ion beam 



sputtering the surface of our nanowire is flat with few nm accuracy [25|. This is quite small 
in absolute terms, particularly keeping in mind that the lattice constant in aluminum is ~ 
0.4 nm. However, approaching the 10 nm scale, a few nm roughness means a lot in relative 
terms. Hence, we cannot exclude the possibility that the specified average effective diameter 
of our nanowires is not very representative: the phase slippage might happen in the weakest 
(thinnest?) parts. To eliminate completely this uncertainty one should study infinitely 
long atomically flat structures, which are quite problematic to fabricate. Nevertheless, the 
main conclusion of the present paper is still valid even assuming the presence of inevitable 
imperfection of the samples. At all realistic parameters characterizing the geometry of the 
possible constrictions, the thermal activation scenario ^ fails to explain the broad R{T) 
dependences observed in sub-15 nm nanowires, while the QPS mechanism ll|] provides 
good agreement with the experiment (Fig. 2). 

In conclusion, we have traced the evolution of the shape of the superconducting transi- 
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tion R{T) in quasi-one-dimensional aluminum nanowires with progressive reduction of their 
cross-sections. For relatively thick samples a^^"^ > 20 nm the shape of the transitions can 
be qualitatively explained in terms of wire inhomogeneity and the model of thermally ac- 
tivated phase slips (TAPS). While in nanowires with effective diameter < 15 nm the R{T) 
dependences are much wider and no reasonable set of fitting parameters can account for 
the TAPS mechanism. The phenomenon is associated with manifestation of quantum fluc- 
tuations of the order parameter. The results are in good agreement with renormalization 



theory [ll|]. The effect of quantum fluctuations should have a universal validity, indicating 
a breakdown of the zero resistance state in quasi-one- dimensional superconductors, setting 
fundamental limitations on miniaturization of nanoelectronic components designed to carry 
a dissipationless supercurrent. 
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